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Abstract 

The aim of this paper is to study the steady states of the mathe- 
matical models with delay kernels which describe pathogen-immune 
dynamics of many kinds of infectious diseases. In the study of mathe- 
matical models of infectious diseases it is important to predict whether 
the infection disappears or the pathogens persist. The delay kernel is 
described by the memory function that reflects the influence of the 
past density of pathogen in the blood and it is given by a nonnega- 
tive bounded function k defined on [0, oo) and is normated. By using 
the coefficient of kernel k, as a bifurcation parameter, the models are 
found to undergo a sequence of Hopf bifurcation. The direction and 
the stability criteria of bifurcation periodic solutions are obtained by 
applying the normal form theory and the center manifold theorems. 
Some numerical simulation examples for justifying the theoretical re- 
sults are also given. 
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1 Introduction 



The purpose of this paper is to study the Hopf bifurcation of pathogen- 
immune dynamics in the steady states of the mathematical model with delay 
kernel . Dynamical systems with delay kernel have been studied for popula- 
tion dynamics and neural networks [2]. 

We introduce a model which describes one of most known infectious dis- 
eases, namely malaria infection. Our model is based on the model from [6]. 
The model, as it was created, without delay, has the feature that the interior 
equilibrium is always asymptotically stable if it exists. For obtaining the 
natural behavior which was observed experimentally, i.e. oscillatory beha- 
vior, it is needed to make some adjustments to the model. One way is to 
introduce some extra terms into the equations for a better description of the 
interaction pathogens-immune system or other way is to introduce the delay. 
Our contribution to the model lies in introduction of the delay kernel, which 
is a natural thing to do, because it is obvious for everyone that biological 
processes do not take place instantaneous. 

System (1) without the last term in the third equation, which implies the 
effect of absorption of the pathogens into uninfected cells, could be used for 
describing another well-known infectious disease, the HIV infection, as was 
done in [3,7,10] or other infectious diseases as hepatitis B virus infection [9] or 
hepatitis C virus infection [7] . Because our model deals with malaria infection 
in what follows we will say some words about this disease. Malaria ranks high 
on the list of world health problems by causing massive human and economic 
loss. The parasite in this disease is called Plasmodium Falciparum and has 
a high virulence which can cause even death. Malaria has some important 
features that should be mentioned: first, the amount of variation in disease 
severity observed in the field is remarkably high, second, the immunity is 
virtually never sufficient to prevent infection and third, transmission intensity 
in the field is highly variable both temporally and geographically [4]. 
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More precisely, consider the following system: 



u 

x(t) — ai — a 2 x(t) — a 3 x(t) J z{s)k{t — s)dt 



— oo 





y(t) = —a4y(t) + a 3 x(t) J z{s)k{t — s)ds (1) 



— oo 





z(t) = a A a^y{t) — a e z(s)k(t — s)ds — a 7 x(t) / z(s)k(t — s)ds, 



where Oj, % = 1, • • • , 7 are positives constants and the delay kernel, k : 
[0, oo) — > [0, oo), k is piecewise continuous, is assumed to satisfy the fol- 
lowing properties: 



oo oo 

J k(s)ds = 1, J sk(s)ds < oo. (2) 



It is also assumed that the system (1) is supplemented with initial con- 
ditions of the form: 

x(0)=x*, y(0) = y*, z(s)=<p 1 (s), se(-oo,0], 
ipi is bounded and continuous on [0, oo). 

The model contains three variables: the density of uninfected cells x, the 

density of infected cells y and the density of pathogens in blood z. Uninfected 

cells are recruited at a constant rate a\ from the source within the body, 

such as the bone marrow and have the natural life expectancy of ^- days. 

Cells are infected by contact with pathogens and turn to infected cells at 
o 

rate a 3 J z(s)k(t — s)ds. Infected cells die at rate a 4 . The death of the cells 

— oo 

results in the release of a$ pathogens per an infected cell and these pathogens 
have a life expectancy of For a 7 = we obtain the model of HIV infection 
with delay kernel and for a 3 = a 7 we obtain the model of malaria infection 
with delay kernel. 
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In system (1), if the delay kernel has the form 

k(s) = 6(s-r), r>0, 

where r is a parameter which denotes the effect of the past memories, then 
system (1) becomes: 

x(t) — ai — a 2 x(t) — a 3 x(t)z(t — r) 

y(t) = -a 4 y(t) + a 3 x(t)z(t - r) 

z(t) = a 4 a 5 y(t) — a e z(t — t) — a 7 x(t)z(t — r) 

x(0)=x*, y(0)=y*, z(s) = Vl (s), s e [— r, 0] 

For r = the system has been studied in [6]. 

In system (1), if the kernel k has the form 

k(s) = qe~ qs , 

called weak kernel, where q is a parameter varying in (0, oo) which denotes 
the decay rate of the effect of the post memories, then system (1) becomes: 

x(t) — ai — a,2x(t) — a 3 x(t)u(t) 
y(t) = -a A y(t) + a 3 x(t)u(t) 

z(t) = a 4 a 5 |/(t) — a 6 u(t) — ajx(t)u(t) (4) 
u{t) = q{z{t) - u{t)) 

x(0)=x*, y(0)=y*, z(s) = z*, u(0) = z* 
In system (1), if the kernel k has the form 

k(s) = q 2 se- qs , 

then system (1) becomes: 

x(t) — ai — Cb2x(t) — a 3 x(t)v(t) 

y (t) = -a A y{t) + a 3 x(t)v(t) 

z(t) = a 4 a 5 y(t) - a 6 v(t) - a 7 x(t)v(t) 

u(t) = q(z(t) - u(t)) { } 

v(t)=q(u(t)-v(t)) 

x(0)=x*, 2/(0) z(0)=z*, u(0) = z*, v(0)=z*. 
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This paper is organized as follows: in Section 2, the local stability pro- 
perty and Hopf bifurcation of models (3), (4), (5) are discussed and some 
sufficient conditions for stability are derived. In Section 3, model (1) con- 
taining the general kernel is further studied and both the direction and the 
stability of Hopf bifurcation are analyzed by the normal form theory and the 
center manifold theorem and some criteria for stability are derived. Then, 
we consider two cases: in the first case k is delta function and in the second 
k = qe~ qs . Numerical simulations will be shown, in order to justify the theo- 
retical results. Finally, conclusions are drawn with further research directions 
given in Section 5. 

2 Local stability analysis and the Hopf bifur- 
cation 

In this section, consider the local stability of the equilibrium solution of 
system (1). From the special nature of the delay kernel (2) embedded in 
system (1), we found out that an equilibrium solution (xq, yo, Zq) of (1) is 
given by the solution of the system: 

ai — a 2 x — a 3 xz = 

a 4U — a 3 xz = (6) 
°4 a 52/ — a 6 z — a 7 xz = 0. 

From (6) it results that if < a 7 < a 3 a 5 , < a 2 a 6 < ai(a 3 a 5 — a 7 ), then 
system (6) has two equilibria. The first one is Xi = 0, 0) and it represents 
the state in which the pathogens are absent. The second is X 2 = (x , y , z ), 
where 

a 6 _ ai(a 3 a 5 - a 7 ) - a 2 G6 _ ai(a 3 a 5 - a 7 ) - a 2 a 6 

— , yo — 7 v , zq — . 

a 3 a5 — a 7 04(0305 — a 7 ) a 3 a6 

The equilibrium X 2 lies in the interior of the first quadrant. Then we say that 
X 2 is an interior equilibrium and represents the state in which the pathogens 
are present. 

In what follows, the equilibrium (xo,yo,^o) is transformed to the origin, so 
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the system (1) becomes: 



u u 
iii(t) — — biUi(t) — b 2 J k(—s)u 3 (t + s)ds — b 3 Ui(t) J k(—s)u 3 (t + s)ds 

— oo — oo 



u 2 (t) = b4Ui(t)-b 5 u 2 (t)+b 2 / k(-s)u 3 (t+s)ds+b 3 ui(t) / k(-s)u 3 (t+s)ds 



o 



-oo 




ii 3 (t) = -b % ui(t) +b 7 u 2 (t) -b 8 J k(-s)u 3 (t+s)ds-b 9 ui(t) j ' k(-s)u 3 (t+s)ds, 



where 



bi = a 2 + a 3 z , b 2 = a 3 x b 3 = a 3 , 6 4 = a 3 z , b 5 = a 4 , 
&6 = 07^0? b 7 = 0405, b$ = + a 7 Xo, bg = a 7 



and 



u^t) = x(t) - x , u 2 (t) = y(t) - y , u 3 (t) = z(t) - z . 
Rewrite system (7) in the following matrix form: 





u 



where 



{t) = Lu{t) + J F(s)u(t + s)ds + H(u(t)), 



u(t) = ( Ul (t),u 2 (t),u 3 (t)) T , 



(7) 



-61 \ / -6 2 

h -b 5 ,F(s) = k(-s) b 2 
-b e b 7 / \ -b 8 

( \ 

— b 3 ui(t) J u 3 (t + s)k(—s)ds 



H(u(t)) 



b 3 u\{t) j u 3 (t + s)k(—s)ds 

—00 


—bgUi(t) f u 3 (t + s)k(—s)ds 

—00 
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(8) 



The associated characteristic equation of the linearized system is given 

by: 

o 

A 3 +p2A 2 +piA + (r 2 A 2 + riA + r ) J k(-s)e Xs ds = 0, (9) 

— oo 

where 

p 2 = b 1 + b 5 , pi = bib 5 , r 2 = b 8 , n = 6 8 (&i + h) ~ - b 2 b 7 
r = b 4 b 7 - b\b 2 b-j - bib 5 b 8 - b 2 b 5 b 6 . 

Proposition 2.1. If k(s) = 5(s — t), then 

(i) The characteristic equation (9) is given by 

A 3 + p 2 X 2 + Pl X + (r 2 A 2 + ri X + r )e- Xr = 0; (10) 

(ii) For t = the characteristic equation (10) is given by 

A 3 + m 2 A 2 + miA + m = 0, 

where 

m 2 =p 2 + r 2 , m 1 =p 1 +r 1 , m = r ; 
(Hi) If r = 0, the equilibrium X 2 is asymptotically stable if and only if 
m 2 > 0, mi > 0, m > 0, m\m 2 — m > 0; 

(iv) For t = tq, given by 

1 rip 2 ^o ~ (^o ~ ^oPi)(?"o - r 2 uJo) 

Tq — clXCtcin ■ 



where uo is the positive root of the equation 

x 6 + nix 4 + n 2 x 2 + ri3 = 0, 

with 

m=p 2 2 - 2p 1 - r 2 , n 2 = p\ - r\ + 2r r 2 , n 3 = -r 2 , 
i/iere is a Hopf bifurcation. 



Proposition 2.2. If k(s) = qe qs , s > 0, q > 0, then 
(i) The characteristic equation (9) is given by: 

A 4 + (p2 + g)A 3 + (pi + q(p 2 + r 2 ))A 2 + g(pi + n)A + r g = 0; 

^ TTie equilibrium X 2 is asymptotically stable if and only if 



(Hi) If there exists q > so that D 3 (q ) = and \ q=qo ^ 0, then a 

Hopf bifurcation occurs at X 2 as q passes through q . 

3 Stability of the bifurcating periodic solu- 
tions: the general kernel case 

In this section, the stability of the bifurcating periodic solutions of system 
(1) with the kernel satisfying (2) is studied. For convenience, in the study of 
the Hopf bifurcation problem, first we transform system (7) into an operator 
equation of the form 



Aj(<7) = ((pi + n)(p 2 + r 2 ) - r )q 2 + ((pi + ri)(p 2 (P2+ 
+ r 2 ) - ri) - 2p 2 r )q + p 2 (pi - r ) > 0; 



ii t = A(fi)u t + Ru 



where 



u = (ui,u 2 ,u 3 ) T , u t = u{t + 6) 



9 G (— oo, 0), fi = a — a 



and operators A and R are defined as 





where 



o 




— oo 
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with L, F defined in (8). 

As in [1], [5], the bifurcating periodic solutions x(t,fi) of (1) are indexed by 
a small parameter e > 0. A solution x(t,fi(e)) has amplitude 0(e), period 
T(e) and nonzero Floquet exponent (5(e) with (3(0) = 0. Under the present 
assumptions, /i, T and (3 have expansions: 



/i = jj, 2 e 2 + /i 4 e 4 + • • • 

T = —(l + T 2 e 2 + T 4 e A + ---) 

(3 = (3 2 e 2 + (3 A e 4 + • • • . 

The sign of \i 2 determines the direction of bifurcation, while (3 2 determines 
the stability of x(t, fi(e)): asymptotically orbitally stable if (3 2 < 0, but 
unstable if (3 2 > 0. 

Next, the question of how to derive the coefficients in these expansions 
is addressed. For the applications from this paper, only fi 2 , r 2 and (3 2 are 
computed here. 

We define the adjoint operator A* of A as: 



A*iP(s) 



- d ^, S G(0,oo) 
L T ij(0) + f F T (s)ij(-s)ds, s = 0, 



where L T and F T are transposes of matrices L and F respectively. 

Note that the operator A depends on the bifurcation parameter a. Ac- 
cording to Propositions 2.1, 2.2, Hopf bifurcation occurs when a passes 
through a . Let /i = a — a . Then, Hopf bifurcation occurs when jj, = 0. It is 
therefore reasonable to assume that (p,v( : [0, oo) — > C 3 . Define the bilinear 
form: 

o e 

< 0, ^ >= W) T m - J J ^ T (£ - e)F(e)<Kt)d£d0. 

6»=-oo^=0 

To determine the Poincare normal form of operator A, we need to calculate 
the eigenvector <p of A associated with eigenvalue Ai = iu and the eigenvec- 
tor <p* of A* associated with eigenvalue A 2 = Ai. 
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Proposition 3.1. (i) The eigenvector <fi of A associated with eigenvalue 
Ai = iu is given by (f)(9) = ve Xl0 , 9 e (—00, 0], where v = (vi,v 2 ,Vs) T 
and 

vi = -b 2 (\i+h)k 1 ,v 2 = & 2 (Ai+&i-& 4 )A;\u3 = (Ai+6i)(Ai+6 5 ), (11) 
where 



k 1 = J k(-s)e XlS ds; 

—00 

(ii) The eigenvector <p* of A* associated with eigenvalue \ 2 = \\ is given 
by 0*(s) = we X2S , s G [0, 00), where w = (wi,w 2 ,W3) T and 



b±b 7 - b 6 (X 2 + 65) 1 A2 + &5 

h{^2 + Oi)r] r] b 7 r] 

b 4 b 7 - 6 6 (A 2 + b 5 )_ _ ,X 2 + b 5 b A b 7 - 6 6 (A 2 + 65) 

v = — m — n-\ — v i + V2 + ( — : (-62 



6 7 (A 2 + 6i) b 7 ' & 7 (A 2 + 6i) 

(12) 



+ b 2 -b 8 ^)k^)v 3 , 

7 



where 



u 

k { ~ 1] = J k(-s)e x ' 2S ds; 



(Hi) We have: 

<0*,0>=1, < > = < >= 0, <0*,0>=1. 

Next, we construct the coordinates of the center of the manifold Q at 
\i = (a = a ) [1], [5]. Let 

z(t) =< <p*,x t > 

w(t,9) =x t -2Re{z(t)<p(9)}. 
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On the center manifold f2 , w(t, 9) = w(z(t), z(t), 9), where 

z 2 z 2 
w(z, z, 9) = w 20 (9) — + w u (9)zz + w 02 (9) — H , 

z and z are the local coordinates of the center manifold Qq in the direction 
of 4> and 0*, respectively. 

For the solution x t G f2 of (1), notice that for /j, — we have: 
z(t) = X 1 z(t)+ < 0*, R(u + 2Re{z(t)(f)(9)}) > 

Rewrite this as 

z(t) = \iz(t) + g(z,z) 

with 

g(z, z) = (j)*(0) T R(w(z, z, 0) + 2Re{z(t)<j>(0)}) 

Further, expand the function g(z,z) on the center manifold fl in powers 
of z and z: 

z 2 ~z 2 z 2 ~z 
g{z,z) = #2oy + 911*2 + 902^ + #2i-^-H 

Proposition 3.2. For the system (1) we have: 

(i) 

g 20 = -263^1^3(^1 -W2~+ ^-w 3 )k l 

03 

gn = -b 3 (viv 3 + v 1 vH)(w 1 -w 2 + -^-w 3 ){k l + k^) (13) 

b 

g 02 = -2b 3 v 1 v^(w 1 -w 2 + — U 3 )A; (_1) ; 

03 



(n) 



(q\ 920 \,0 920 _ \ 2 g J-, 2Ai6» 

w 20 {9) = —ve 1 - -rr-ve + E x e 
\\ 0A1 

tn\ 9ll A i0 9ll— \o8 , 7-1 

wii(0) = — we Aiy - — ve A2V + £ 2 , 
Ai Ai 
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where Ei, E 2 are the solutions of the following system: 

(A + k {2) B - 2X 1 I)E 1 = 63^1^3^(1, -1, -b 9 /b 2 3 ) T 

(A + B)E 2 = b 3 (v lV3 + v.v^k 1 + fc(" 1 ))(l, -1, -h/blf 


k w = J k(-s)e XlS ds,kW = J k(-s)e 2XlS ds; 



(Hi) 



u 

b f 

g 21 = -26 3 (uJi -ufi+ yW 3 )(vi J k(-s)w 311 (s)ds+ 







(14) 



I r 1 

+ -vi / k(-s)w 320 {s)ds + vsWmiOjk 1 + -^3^120 (0)£; (_1) ), 



-00 



with w 20 (9)=(w 120 (9), w 220 (9), w 320 (9)) T and w 11 (9) = ( y w 111 (9), w 211 (9) 
w 311 (9)) T - 

Therefore, we will compute the following parameters: 



2 1| |2\ , 921 



ci(0) = ^(^20^11 - 2 |^ii I - -\g 02 1 ) + — 
Re ci(0) 



i?e A'(ao) 

Jmci(O) + /i 2 /m A'(a ) 
j 2 _ 

^0 

p 2 = 2Re Cl (0) 

T = — (1 + T 2 e 2 + 0(5 4 )), e 2 = + 0(a - a ) 2 . 

We have: 

Theorem 3.1. The sign of fx 2 determines the directions of the Hopf bi- 
furcations: if /jl 2 > 0(< 0) the Hopf bifurcation is supercritical (subcritical) 
and the bifurcating periodic solutions exist for a > a (< a ). The sign of f3 2 
determines the stability of the bifurcation periodic solutions. They are both 
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asymptotically orbitally stable if (3 2 < 0, but unstable if (3 2 > 0. T 2 deter- 
mines the period of the bifurcating periodic solutions: the period increases 
(decreases) if T 2 > 0(< 0). 

If k(s) = 5(s - r),r > 0, then k 1 = e AlT , fc(" 1 ) = e AaT , k 2 = e 2AlT and 
a = t , where r is given by (20). In this case 



t=t , A=Ai — 

Ai(r 2 A? + nAi + r ) 



(3A? + 2p 2 Ai + pi)e A i r o + 2r 2 A x + n - (r 2 A? + nAi + r )r 

If fc(s) = ge-^,s > 0,g > 0, then A; 1 = k^ = ^and a = q , 

where go satisfies D 3 (q ) = (D 3 (qo) from Proposition 2.2), Ai = iuo, A 2 = Ai 
and u>o is given by 

2 9o(pi + ri) 

w o ; • 

P2 + qo 

In this case 

A'(<?o) = -J-\q=q ,X=X 1 = 



dq 



A? + (pa + r 2 )A? + (p! + n)Ai + r 



4A? + 3(p 2 + q )\j + 2(p! + q (p 2 + r 2 ))Ai + g (Pi + n) ' 

From (11), (12), (13), (14) it results: 

Proposition 3.3. If k(s) = qe~ qs , s > 0, q > 0, then for system (1) we 
have: 

b 2 (Xi + b 5 )q b 2 (Xi + h - b 4 )q 

vi = t — ■ , v 2 = — ■ , v 3 = (Ai + 6i)(Ai + 6 5 ), 

Ai + g Ai + go 



6 4 6 7 - 6 6 (A 2 + 65) 1 A 2 + b 5 

Wl = T~T\ 71—^ > w 2 = ~, ^3 



b 7 (X 2 + b 1 )r] 77' b 7 r] 



b 4 b T - b 6 (X 2 + 65) _ ,A 2 + fc 5 b 4 b 7 - 6 6 (A 2 + 65) 

" = 6 7 (A 2 + 61) + " 2 + ( ~ MA^TV) 

. , , A 2 + 6 5 g 

+ 6 2 - 6 8 — )- — ■ — )v 3 

b 7 A 2 + g 
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#20 = -2b 3 V 1 V 3 (lU 1 -W 2 + 7-W 3 ) ; 

03 Ai + g 

, /_ . — w_ _ . b 9 _ q g 

9n = -h{viv 3 + ViV 3 ){w 1 -w 2 + t-w 3 ){- — ; h t — ; ) 

t> 3 Ai + qo a 2 + go 

, /_ _ . h— x go 

g 02 = -263^1^3(^1 - w 2 + r w 3 jT — ; 

03 A 2 + go 



5-21 = -26 3 - w 2 + t-w 3 ) «i(t7^ — i - VT^ — ; \ v 3 + £32) 

03 Ai(Ai + go) Ai(A 2 + g ) 

, !- / fl'aogo # 2 o<?o — . 7-1 go x 

2 Ai(Ai + g ) 3Ai(A 2 + g ) 2Ai + g 

. f 9u 9n- . P x go 1- f 9wvi 9 20 Vi , p g2 ^ 

where E x = (E n , E 21 , E 31 ) T , E 2 = (E 12 , E 22 , E 32 ) T satisfy the systems: 
^ + ^—B - 2X 1 I)E 1 = b 3Vl v 3 -^—(l, -1, -|f 



2Ai + g 2Ai + g b 



(A + B)E 2 = b 3 (v lV3 + Ul u 3 )g (—L_ + — |— )(1, -1, -J|) T . 

Ai + go a 2 + go o 3 



4 Numerical simulations 

For numerical simulations, we use Maple 9.5. We consider system (3) with 
a,\ = 2, a 2 = 0.02, 03 = 0.5, 04 = 2, 05 = 1.5, a§ = 0,03, a 7 = 0.5. The 
equilibrium point is: Xo = 0.12, y = 0.99, Zo = 33,29. 

In the first case k(s) = 5(s — r), we obtain: r = 0.8975032747, w = 
1.140149275, #20 = -0.1032034629 + 0.8542910112i, g u = 0.6157126005 - 
0.5196432315^,^2 = -0.1015556589+0.8738154402i, g 21 = -0.2181868792+ 
0.9831278086i, ci(0) = -0.1197163151 + 0.2827563802i, // 2 = 0.2857035314, 
P 2 = -0.2394326302, T 2 = 0.1030243826. 

Because fj, 2 > 0, the Hopf bifurcation is supercritical for r > r ; as 
(5 2 < the bifurcating periodic solution is asymptotically orbitally stable; as 
T 2 > the period increases. Here are the computer simulations. Fig.l rep- 
resents time vs uninfected cells (t,x(t)), Fig.2 time vs infected cells (t,y(t)), 



14 



Fig. 3 time vs pathogens in blood (t,z(t)), Fig. 4 pathogens vs infected cells 
(z(t),x(t)), Fig.5 pathogens vs infected cells (z(t),y(t)), Fig.6 uninfected cells 
vs infected cells (x(t),y(t)). 
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Fig.5.Phaseplot(z(t),y(t)) 



Fig.6.Phaseplot(x(t),y(t)) 




In the second case k(s) = qe qs , we obtain: q = 0.1881852832, uj = 
0.1872904846, g 20 = -0.1146541958 - 0.3262916164i, g u = 0.2100798781 - 
0. 25546255 Ui, g 02 = -0.1034928343+0.5915640328z, g 21 = -0.1382064847- 
0.2196567760i, ci(0) = -0.6909931661-0.1103278787i, fi 2 = -0.4372331513, 
(3 2 = -0.1381986332, T 2 = 0.8650529366. 

Because fi 2 < 0, the Hopf bifurcation is sub critical for q > q ; as [3 2 < 
the bifurcating periodic solution is asymptotically orbitally stable; as T 2 > 
the period increase. Here are the computer simulations: Fig. 7 represents time 
vs uninfected cells (t,x(t)), Fig.8 time vs infected cells (t,y(t)), Fig. 9 time 
vs pathogens in blood (t, z(t)), Fig.10 pathogens vs infected cells (z(t),x(t)), 
Fig. 11 pathogens vs infected cells (z(t),y(t)), Fig. 12 uninfected cells vs in- 
fected cells (x(t),y(t)). 



Fig.7.(t,x(t» Fig.8. (t,y(t» 




20 40 60 80 20 40 60 80 

t t 
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Fig.9.(t,z(t)) 



Fig.10.Phaseplot(x(t),z(t)) 




32.400002 32.799999 33.200001 33.599998 34 0.1297 0.1298 0.1299 0.13 0.1301 0.1302 0.1303 



5 Conclusion 

In this paper we introduce a model which describes infectious diseases and 
malaria infection with delay kernel. By using the average time delay as a 
parameter, it has been proved that the Hopf bifurcation occurs when this 
parameter passes through a critical value. A similarly study will be made for 
the models which describe infectious diseases, who takes into account immune 
response against pathogens and the effect of involvement. This study will be 
made in a future work. 
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